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> '■ Abstract 

We show that the interplay between the LHC and the e + e~ International Linear 
Collider (ILC) with yfs = 500 GeV might be crucial for the discrimination between 
the minimal and next-to-minimal supersymmetric standard model. We present an 
NMSSM scenario that cannot be distinguished from the MSSM by cross sections 
and mass measurements if only the light neutralinos and the lightest chargino are 
kinematically accessible, even if one of the neutralinos has a significant singlino 
component. Mass predictions for the heavier neutralinos from the ILC analysis 
and their observation at the LHC lead to an identification of the neutralino mix- 
ing character and the underlying supersymmetric model in a combined LHC /ILC 
analysis. In our numerical example we include errors in the mass measurements 
and use standard methods of supersymmetric parameter determination. 



1 Introduction 



The International Linear Collider (ILC) is intended to start with an energy of ^fs = 
500 GeV, which will be upgraded to about 1 TeV Q. The physics program includes the 
complement of the results at the Large Hadron Collider (LHC) |2] and the discovery of 
possible signatures of physics beyond the standard model (SM). Since supersymmetry 
is often considered as the main candidate for new physics, the precise resolution of the 
supersymmetric model is one of the essential goals of the ILC 0. 

An interesting possibility for the determination of the supersymmetric model is to 
study the gaugino/higgsino particles, which are expected to be among the lightest su- 
persymmetric particles. In this paper we consider two basic supersymmetric models: 
the minimal supersymmetric standard model (MSSM) and the next-to-minimal super- 
symmetric standard model (NMSSM). The MSSM contains four neutralinos Xi/ which 
are mixtures of the photino, zino and neutral higgsinos, and two charginos xf, be- 
ing mixtures of wino and charged higgsino. The neutralino/chargino sector depends 
at tree level on four parameters: the U(l) and SU(2) gaugino masses M\ and M 2 , the 
higgsino mass parameter n, and the ratio tan (3 of the vacuum expectation values of 
the Higgs fields. For the determination of these parameters, straightforward strategies 
IIHEI have been worked out even if only the light neutralinos and charginos Xi> X2 an d 
xf are kinematically accessible at the first stage of the ILC [6|. 

The NMSSM |7[ is the simplest extension of the MSSM by an additional Higgs sin- 
glet field. New parameters in the neutralino sector are the vacuum expectation value 
x of the singlet field and the trilinear couplings A and k in the superpotential, where 
the product Ax = /x eff replaces the /^-parameter of the MSSM ||HH3. The additional fifth 
neutralino may significantly change the phenomenology of the neutralino sector. In 
scenarios where the lightest supersymmetric particle is a nearly pure singlino, the ex- 
istence of displaced vertices leads to a particularly interesting experimental signature 
ITOllllll . If the complete neutralino sector is accessible, the sum rules for the produc- 
tion cross sections show a different energy dependence in MSSM and NMSSM |6[. The 
NMSSM Higgs sector also depends on the soft scalar mass parameters A\ and A K . It 
contains five physical neutral Higgs bosons, three scalars Sj and two pseudoscalars 
Pj. In the following we will assume that only the light neutralinos, charginos and the 
lightest Higgs boson Si are accessible at the first stage of the linear collider. 

It has already been worked out that there exist MSSM and NMSSM scenarios with 
the same mass spectra of the light neutralinos but different neutralino mixing. In this 
case beam polarization is crucial for distinguishing the two models [12J. We present 
a scenario where all kinematically accessible neutralinos and charginos have simi- 
lar masses and almost identical cross sections, within experimental errors, in MSSM 
and NMSSM. Although the second lightest neutralino in the NMSSM has a significant 
singlino component, the models cannot be distinguished by precision measurements 
of masses and cross sections at the LHC or at the ILC with ^fs = 500 GeV. While the 
branching ratios for the neutralino decays into Higgs bosons may give first evidence, 
final identification of the underlying model requires precision measurements of the 
masses and mixing characters of the heavier neutralinos by combined analyses of LHC 
and ILC500 results and by exploiting the ILC 65 ^ ' option with y/s = 650 GeV at the cost 
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of reduced luminosity. The ILC650 option reflects the high flexibility of the ILC and 
does not require any specific setups in the experiment. Of course, the second energy 
stage of the ILC with energies up to 1 TeV also allows a discrimination. Here, however, 
we focus on the potential of the ILC 50 o/65o- 

This paper is organized as follows: in Section 2 we describe in detail our strategy for 
the discrimination between NMSSM and MSSM. Our analysis is based on the scenario 
presented in Section 3. The following sections 4, 5 and 6 contain the numerical results 
for the parameter determination at ILC 50 o/ LHC and ILC 650 . We estimate errors in the 
mass measurements and include statistical errors for the cross sections. Fit programs, 
which exist so far only for the MSSM |13|, are not used. Special attention is given to 
the interplay in combined analyses of the experimental results at ILC and LHC. 

2 Strategy 

We consider an NMSSM scenario where in the gaugino/higgsino sector only the 
masses of the light chargino xt an d the light neutralinos Xi an d X2 are accessible at 
the ILC500 We calculate the masses of the charginos and neutralinos and the cross 
sections for the pair production of the light chargino e + e~ — > XiXi an d for the associ- 
ated production of the light neutralinos e + e — > X1X2 with polarized and unpolarized 
beams. 

The masses and cross sections provide the experimental input for deriving the su- 
persymmetric parameters within the MSSM using standard methods 

• We assume an uncertainty of 1-2% for the masses m^±, m^o, m^o, mp e , and 

m e - H . In the cross sections for e + e~ — > XiXi an d e + e~ — > X1X2 we include a 
statistical error of one standard deviation and take into account the error due to 
the polarization uncertainty of AV e ± /V e ± = 0.5%. 

• From the chargino mass m^± and the cross section e + e~ — > XiXi measured at two 
energies, yjs = 400 GeV and 500 GeV, we determine bounds for the elements U\\ 
and V\\ of the chargino mixing matrices. Polarized beams allow the resolution of 
ambiguities and the improvement of the accuracy. 

• Using the mixing matrix elements U\\ and V\\, the masses rn^±, m^o and m^o, and 
the cross sections for e + e~ — > xixtf we derive constraints for the parameters M t , 

M 2/ \i and tan (3. 

• After the determination of the fundamental MSSM parameters we calculate the 
heavy chargino and neutralino masses and compare them with LHC analyses 

m 

The experimental results from the ILC with y/s = 400 GeV and 500 GeV lead to consis- 
tent MSSM parameters despite the fact that we start with an NMSSM scenario where 
the light neutralino x° nas a significant singlino component. In the considered scenario 
the ILC500 does not allow a discrimination between MSSM and NMSSM by mass and 
cross section measurements. Then the identification of the underlying supersymmetric 
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model is performed in combined analyses at the LHC or the ILCg 50 1,/3 by the measure- 
ment of the masses of the heavy neutralinos and charginos. 

3 Scenario 

We discuss an NMSSM scenario where the lightest visible neutralino x° is accessible 
at the ILC500 and has a significant singlino component. The supersymmetric param- 
eters and the resulting masses and neutralino mixing states are given in Table [TJ For 
the basis in the neutralino system we use the conventions of IT5lll6l . The hierarchy 
Mi > M 2 of the U(l) and SU(2) mass parameters causes an approximate mass degener- 
ation of the lightest neutralino Xi, which is assumed to be the lightest supersymmetric 
particle (LSP), and the light chargino Xi that is typical for a minimal anomaly medi- 
ated SUSY-breaking mechanism (mAMSB) |17[. Even small mass differences may be 
resolved experimentally with good accuracy at the LHC [ 18 [ and by applying the ISR 
method |19[ at the linear collider |20|. One should note, however, that we do not use 
any assumptions on possible SUSY-breaking mechanisms for the following parame- 
ter determination. For comparison Table [T] also contains the parameters of the MSSM 
scenario derived in Section 4, which leads to indistinguishable mass spectra and cross 
sections within the experimental errors at the ILC 50 o- All parameters are in agreement 
with bounds from dark matter constraints 1211 . 

Table|2]shows the masses and mixings of the NMSSM Higgs bosons for parameters 
A\ and A K within the theoretical and experimental constraints in the NMSSM scenario 

2740 GeV < A x < 5465 GeV, (1) 
-553 GeV < A K < GeV. (2) 

Here the light scalar Si has MSSM-like character while the second lightest scalar S2 and 
the light pseudoscalar Pi are almost pure singlets with masses larger than 200 GeV for 
—443 GeV < A K < —91 GeV. Within this A K region direct Higgs production does not 
allow the identification of the NMSSM [22 1 since scalar and pseudoscalar Higgs bosons 
with dominant singlet character escape detection. Detailed studies on the detection of 
light singlet dominated Higgs bosons at the LHC and ILC can be found in |23|, and of 
light pseudoscalars in scenarios with dominant decay S\ — > Pi Pi at LHC, ILC and the 
photon collider in |24[. 

Due to the mass difference between the light neutralinos in our scenario the decays 
of X2 an d X3 u^to S2 and Pi are kinematically forbidden. The branching ratio for the 
decay of X2 too the lightest Higgs boson Si is enhanced by a factor two to about 20% in 
the NMSSM scenario compared to the MSSM. Therefore precise measurements of the 
neutralino decays into Higgs bosons could provide first evidence for the underlying 
model that has to be confirmed by the combined LHC and ILC analyses described in 
the next sections. 
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NMSSM 



MSSM 



Mi 
M 2 
tan (3 

A 

x 

K 



360 GeV 
147 GeV 
10 

0.5 
915 GeV 
0.2 



375 GeV 
152 GeV 
8 

360 GeV 



A2 
A3 

m xl 
m x% 



138 GeV 
337 GeV 
367 GeV 
468 GeV 
499 GeV 



138 GeV 
344 GeV 
366 GeV 
410 GeV 



m xt 
m xt 



139 GeV 
474 GeV 



139 GeV 
383 GeV 



m 



''L 



m 



'■R 



240 GeV 
220 GeV 
226 GeV 



240 GeV 
220 GeV 
226 GeV 



Xi mixing state 
X2 mixing state 
xl mixing state 
xl mixing state 
xl mixing state 



(-0.02, 0.97, -0.20, 0.09, -0.07) 
(0.62, 0.14, 0.25, -0.31, 0.65) 

(-0.75, 0.04, 0.006, -0.12, 0.65) 
(-0.03, 0.08, 0.70, 0.70, 0.08) 
(0.21, -0.16, -0.64, 0.62, 0.37) 



(0.03, -0.96, 0.26, -0.13) 
(0.72, 0.22, 0.48, -0.46) 
(-0.04, 0.10, -0.70, -0.71) 
(-0.70, 0.18, 0.47, -0.52) 



Table 1: Parameters, mass eigenvalues and neutralino mixing states in the basis 
(B°, W°, Hi,H$,S) in the NMSSM scenario. For comparison the table also contains 
an MSSM scenario (neutralino basis (B°, W°, H®, H%)) with a similar mass spectrum of 
the light neutralinos and charginos but different neutralino mixing. 
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A x 


4000 GeV 


A K 


-200 GeV 


m Sl 


124 GeV 


m S2 


311 GeV 


m S3 


4347 GeV 


m Pl 


335 GeV 


mp 2 


4347 GeV 


m H ± 


4346 GeV 


S\ mixing state 


(-0.9947, -0.0990, 0.0292) 


S 2 mixing state 


(-0.0271, -0.0222, -0.9994) 


53 mixing state 


(0.0996, -0.9948, 0.0193) 


Pi mixing state 


(0.0016, 0.0164, -0.9999) 


P 2 mixing state 


(0.0995, 0.9949, 0.0165) 



Table 2: The parameters A\ and A K of the NMSSM Higgs sector and the resulting 
masses and mixing states of the Higgs particles according to NMHDECAY II15II for the 
parameters in the neutralino sector as given in Tableland Mq = My = M D = 800 GeV 

and A t = A b = 1500 GeV. 

4 Parameter determination at the Linear Collider 

In order to analyse the possible parameter reconstruction at the ILC 50 o, we start with 
the NMSSM scenario of Table [T] and calculate the resulting cross sections for neu- 
tralino and chargino production for unpolarized beams and beam polarizations of 
(p e _,p e+ ) = (-90%, +60%) and (+90%, -60%) and assume AV e -/V e - = AV e +/P e+ = 
0.5% |25|. At y/s = 400 GeV chargino pairs xtxT are produced, at y/s — 500 GeV the 
pairs X1X1 an d X1X2 are visible. Now the well-known strategies for determining the 
fundamental MSSM parameters in the gaugino/higgsino sector |6] are applied, tak- 
ing into account experimental errors of 1.5% (2%) in the mass measurements of X23/ 
£l,r, v (Xir Xi) an d one standard deviation statistical errors for the cross sections 1141 
as given in Table |3j The statistical errors are based on a total luminosity of 500 fb _1 , 
which results in J C = 100 fb^ 1 for each of the possible five polarization configurations 
(Ve-,V e +) = (0,0), (-90%, +60%), (-90%, -60%), (+90%, +60%) and (+90%, -60%). 

We determine the components of the chargino eigenstates Uu and Vn |26| from 
the measured polarized cross sections, which are bilinear quadratic functions of U\ x 
and V^i and depend also on m -± and . The expressions for the chargino cross sec- 
tions are analytically inverted jl4l . In order to resolve ambiguities we assume that 
the process has been measured at two different energies, y/s = 400 and 500 GeV, and 
determine the components of the chargino eigenstates, which are consistent with the 
measured cross sections within the estimated error bounds of Table 

Uli = [0.84, 1.0], V11 = [0.83, 1.0]. (3) 
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m^o 1 GeV 
m^o 1 GeV 
m -± / GeV 

A1 


138 + 2.8 
337 + 5.1 

139 + 2.8 


rrie L 1 GeV 
rrie R 1 GeV 
mo J GeV 


240 ± 3.6 
220 + 3.3 
226 ± 3.4 


(r(e + e- -> xfx?)/fo at = 400 GeV 


Unpolarized beams 

(V e -,V e +) = (-90%, +60%) 
(V e -,V e+ ) = (+90%, -60%) 


323.9 + 33.5 
984.0 + 101.6 
13.6 + 1.6 


a(e + e- -> xfx^/fb at = 500 GeV 


Unpolarized beams 

(V e -,V e+ ) = (-90%, +60%) 
(V e -,V e +) = (+90%, -60%) 


287.5 + 16.5 
873.9 + 50.1 
11.7+1.2 


a(e + e- -> x?x§)/fb at = 500 GeV 


Unpolarized beams 

(V e -,V e +) = (-90%, +60%) 
(V e -,V e+ ) = (+90%, -60%) 


4.0 + 1.2 
12.1 + 3.8 
0.2 + 0.1 



Table 3: Masses with 1.5% (X23, &l,r, v) and 2% (x\, Xx) uncertainty and cross sections 
with an error composed of the error due to the mass and polarization uncertainties and 
one standard deviation statistical error based on J £ = 100 fb 1 , for both unpolarized 
beams and polarized beams with {V e -,V e +) = (+90%, ±60%) and AV e ±/V e ± = 0.5%, 
in analogy to the study in [14J. 



We then exploit additionally the polarized neutralino cross sections a(e + e~ — > X1X2) a * 
y/s = 500 GeV and the measured masses m^o, rn^o, mg L and vrie R in order to determine 
the parameters M x , M 2 , /x and tan (3: 

M x = 377 ±42 GeV, (4) 

M 2 = 150 ±20 GeV, (5) 

H = 450 ±100 GeV, (6) 

tan/3 = [1,30]. (7) 



Note that, in our mAMSB-inspired scenario with M\ > M 2 , the crucial observable to 
determine the parameter Mi is rn^o, which can be clearly seen in Fig. [l] [27J. Since the 
heavier neutralino and chargino states are not produced, the parameters /i and tan (5 
are determined with a considerable uncertainty 

In conclusion the measured masses and cross sections as experimental input from 
the ILC 50 o are compatible with the MSSM. For the allowed ranges ®-0 of the MSSM 
parameters we now predict the heavier neutralino and chargino masses within the 
MSSM: 

m^o = 443 ±107 GeV, (8) 
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Figure 1: The neutralino masses m^o, % = 1, . . . , 5 as a function of Mi in the NMSSM for 
the scenario of Table [TJ 

rn^o = 490±110GeV, (9) 
m-± = 475±125GeV. (10) 

We emphasize that although we started with an NMSSM scenario where x° an d X° 
have large singlino admixtures (see TableHJ), the strategy [5, 6 J to determine the MSSM 
parameters is successful. All predictions for the heavier gaugino/higgsino masses are 
consistent in both models. In order to distinguish between the models, the measure- 
ment of the heavier neutralino and chargino masses at the LHC or with the ILC 6 5o 
option is necessary. In our NMSSM scenario the gaugino admixture of probably al- 
lows its mass measurement at the LHC, whereas in the MSSM the neutralino xt always 
has a strong higgsino character within the predicted parameter ranges, eqs. ©-0/ and 
would hardly be visible at the LHC. 



5 SUSY searches at the LHC 

The LHC has an excellent potential to discover supersymmetric particles. The expected 
large cross sections for squark and gluino production at the LHC give access to a large 
spectrum of coloured as well as non-coloured supersymmetric particles via the cascade 
decays. Simulations |28[ exist e.g. for the mSUGRA point SPSla 1291 . It has also been 
worked out how to measure the masses of the heavier gauginos in cascade decays at 
the LHC |30J. Mass predictions from the ILC analysis lead to an increase of statistical 
sensitivity for the LHC analysis and open the possibility of identifying heavier gaugi- 
nos even in marginal signals in the squark cascades [14] . However, since higgsino-like 
charginos and neutralinos do not couple to squarks, their detection via cascade decays 
is not possible. Direct pair production of charginos and neutralinos that occurs via 
Drell-Yan production is strongly suppressed for masses larger than about 300 GeV as 
well as single neutralino production [31 J. 

The mAMSB-inspired scenario of Table [l] contains nearly mass degenerate Xi an d 
Xi with an LSP-NLSP mass difference in the difficult region for hadron colliders of 
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200 MeV < m~± — m^o < 2 GeV because the short-lived xf is n °t directly observable 
and its SM decay particles are softly emitted and suffer from a large background l!T8ll . 
Applying specific cuts and using suitable kinematical variables, however, leads to a 
rather accurate measurement of the mass difference. A gain in accuracy can be ob- 
tained using the LSP and NLSP mass from ILC measurements. Since there exists no 
simulation for our chosen parameter point, we assume a precision for the heavier neu- 
tralino masses as in |30|, motivated by supersymmetric searches in AMSB scenarios at 
the LHC |18|. A more detailed analysis will be performed in a forthcoming work II32II . 

The neutralinos x_2 an d X3 have a large bino-admixture and therefore appear in 
the squark decay cascades. The dominant decay mode of x\ nas a branching ratio 
BR(x° 2 -> xfW T ) ~ 50%, while for the xl decays BR(x° 3 -> £% R i T ) ~ 45% is largest. 
Since the heavier neutralinos, xlr X5/ are mainly higgsino-like, no visible edges from 
these particles occur in the cascades. It is expected to see the edges for X2 ~ > £r^ t > 
xl^ Xl - a nd for xl - If^ E3E|. 

With a precise mass measurement of XirXv ^l,r a nd v from the ILC 50 o analysis, a 
clear identification and separation of the edges of the two gauginos at the LHC is pos- 
sible without imposing specific model assumptions. We therefore assume a precision 
of about 2% for the measurement of m^o, in analogy to 1 30 ] : 

rn^o = 367 ± 7 GeV. (11) 

Here the error assumption is rather conservative since contrary to the scenarios in [ 30 1 
the branching ratio BR(q R — > xll) x BR(x\ — > if R t T ) is very large |33). 

The precise mass measurement of Xz is compatible with the mass predictions of the 
ILC500 However, it is still to confirm that the measured particle is indeed the X3 an d 
not the xl- The predicted X3 i n the MSSM is nearly a pure higgsino, which is typical for 
the constrained MSSM, and does not couple in the cascade decays, while the xi with 
a sufficiently large gaugino component could be measured in the cascades I14U30I . In 
order to identify the particle and finally determine the masses and mixing characters 
of the heavy neutralinos one has to discuss the following cases: 

• Assuming the measured particle to be X3 an d feeding it back in the ILC analysis 
lead to improved parameter determination and mass predictions for m^o and 
m-±. Using eq. ((TT) for the ILC 50 o analysis leads in our case, after rechecking 
with the allowed cross sections of X1X2 an d X1X1 production, to the precise mass 
predictions: 

m^o = [384, 393] GeV and m-± = [360, 380] GeV. (12) 

• Assuming the measured particle to be x°i an d feeding it back in the parameter 
determination of the ILC analysis lead to inconsistency with the measured cross 
sections of X1X2 an d X1X1 production. 

The combined LHC/ILC 50 o analysis therefore leads to a consistent interpretation 
of the measured particles in the cascades. However, a neutralino X3 with sufficiently 
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large gaugino admixture to couple to squarks is incompatible with the allowed pa- 
rameter ranges of eqs. ©-(Hi in the MSSM. One should note that this result has been 
derived within the chargino/neutralino sector of the CP-conserving unconstrained 
MSSM, without any assumption concerning the SUSY-breaking mechanism. 

We point out that a measurement of the neutralino masses m^a, rn^o, m^o which 
can take place at the LHC alone is not sufficient to distinguish the SUSY models since 
rather similar mass spectra may exist. Both a model-independent analysis using cross 
section and masses at the ILC and the precisely measured masses from the cascade 
analysis at the LHC are needed for the identification of the particles. Furthermore, 
the results from the LHC and the ILC 50 o analyses and the precise predictions for the 
missing chargino/ neutralino masses, eq. ((12)) , motivates to apply immediately the low- 
luminosity but higher-energy option ILC^ 1 ^ 3 , which finally could lead to the correct 
identification of the underlying model. 



6 The ILC^ 1 3 option 

The ILC design provides the option to enhance instantaneously the energy at the cost 
of reduced luminosity. A centre-of-mass energy of ^fs = 650 GeV at a third of the 
luminosity, i.e. J £ = 33 fb 1 per polarization configuration, does not require hardware 
changes in the experimental set-up but only running time. 

The inconsistency of the MSSM analyses, the predicted higgsino-like x® an d the 
measurement of this neutralino at the LHC, described in the previous sections may 
motivate this higher-energy option. The expected polarized and unpolarized cross 
sections, including the statistical error on the basis of one third of the luminosity of 
the ILC500/ are given in Table HI The neutralino x° as well as the higgsino-like heavy 
neutralino X4 an d the chargino xf are now accessible at the ILCg 5 ^ 1//3 . The production 
of X1X3 an d X1X4 leads to promising rates at the 650 GeV option while the heaviest 
neutralino X5 is still too close to the threshold and the production cross sections are 
probably too small to be observed. The precise mass measurements of x\ an d xt l ea d 
to a clear identification of the supersymmetric model by comparing them with the 
predictions of eq. (H2J. The accuracy of fj, e s = Xx is improved by the measurements of 
the cross sections for the associated X1X2 production shown in Tabled 

Finally the masses m^o, m^o, m^o, m^o and rn^±, m^±, as we ^ as * ne cr °ss sections, 
constitute the observables for a fit of the parameters Mi, M 2 , tan (3, X and /i c ff = Xx, kx 
of the gaugino /higgsino sector, which is, however, beyond the scope of this paper. Our 
analysis shows that the interplay between the two experiments is crucial for the deter- 
mination of the supersymmetric model. It motivates the use of the low-luminosity, 
y/s = 650 GeV, option of the ILC in order to resolve model ambiguities even at an early 
stage of the experiment and outlines future search strategies at an upgraded LHC as 
well as at the ILC at 1 TeV 
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a{e + e~ -> x?x°)/fb at y/s = 650 GeV 


J = 3 


j = A 


j = 5 


Unpolarized beams 


12.2 + 0.6 


5.5 + 0.4 


< 0.02 


(V e -,V e +) = (-90%, +60%) 


36.9 + 1.1 


14.8 + 0.7 


< 0.07 


(V e -,V e +) = (+90%, -60%) 


0.6 + 0.1 


2.2 + 0.3 


< 0.01 



Table 4: Expected cross sections for the associated production of the heavier neutrali- 
nos in the NMSSM scenario of Table [T]f or the ILCg 5 Q 3 option with one sigma statistical 
error based on £ = 33 fb 1 for both unpolarized and polarized beams. 





a{e + e- -> XiXz)/^ at = 650 GeV 


Unpolarized beams 


2.4 + 0.3 


(V e -,V e +) = (-90%, +60%) 


5.8 + 0.4 


(V e -,V e +) = (+90%, -60%) 


1.6 + 0.2 



Table 5: Expected cross sections for the associated production of the charginos in the 
NMSSM scenario of Table [l] for the ILC^ 1 ^ 3 option with one sigma statistical error 
based on C = 33 fb -1 for both unpolarized and polarized beams. 

7 Conclusion 

We have presented a scenario in the next-to-minimal supersymmetric standard model 
(NMSSM) that cannot be distinguished from the MSSM at the first stage of the Inter- 
national Linear Collider with y/s = 500 GeV by cross section and mass measurements. 
Although a light neutralino has a significant singlino component in the NMSSM, the 
masses of the accessible light neutralinos and charginos, as well as the production cross 
sections, are identical in the two models within experimental errors. The measure- 
ment of the masses of the heavier neutralinos and charginos in combined analyses of 
the experimental results at the LHC and at the higher-energy option ILC 650 with 
y/s = 650 GeV leads to a clear identification of the supersymmetric model. 

Acknowledgements 

We are very grateful to G. Polesello and P. Richardson for constructive discussions. 
We also thank J. Kalinowski for his critical reading of the manuscript. S.H. is sup- 
ported by the Goran Gustafsson Foundation. This work is supported by the European 
Community's Human Potential Programme under contract HPRN-CT-2000-00149, by 
the 'Fonds zur Forderung der wissenschaftlichen Forschung' of Austria, FWF Project 
No. P16592-N02, and by the Deutsche Forschungsgemeinschaft (DFG) under contract 
No. FR 1064/5-2. 



11 



References 

[1] J.-E. Augustin et al., Executive Summary, International Technical Review Panel, 
August 2004, see also http : / / www . interactions . org/linearcollider/. 

[2] J. G. Branson, D. Denegri, I. Hinchliffe, E Gianotti, E E. Paige and P. Sphicas [AT- 
LAS and CMS Collaborations], Eur. Phys. J. direct C 4 (2002) Nl. 

[3] T Abe et al. [American Linear Collider Working Group Collaboration], 
hep-ex/0106056: 

TESLA Technical Design Report, Part III, eds. R.-D. Heuer, D. Miller, E Richard 
and P. Zerwas, hep-ph/0106315: 

K. Abe et al. [ACFA Linear Collider Working Group Collaboration], 
|hep-ph/0109166| 

[4] T Tsukamoto, K. Fujii, H. Murayama, M. Yamaguchi and Y. Okada, Phys. Rev. D 
51 (1995)3153; 

J. L. Feng, M. E. P eskin, H. Murayama and X. Tata, Phys. Rev. D 52 (1995) 1418 
piep-ph /95022601. 

[5] S. Y. Choi, A. Djouadi, M. Guchait, J. Kalinowski, H. S. Song and P. M. Zerwas, 
Eur. Phys. J. C 14 (2000) 535 [hep-ph/ 0002033 1; 

S. Y. Choi, A. Djouadi, H. K. Dreiner, J. Kalinowski and P. M. Zerwas, Eur. Phys. J. 
C 7 (1999) 123 |hep-ph/9806279|. 

[6] S. Y. Choi, J. Kalinowski, G. Moortgat-Pick and P. M. Zerwas, Eur. Phys. J. 
C 22 (2001) 563 [hep-ph/0108117|; Addendum, Eur. Phys. J. C 23 (2002) 769 
piep-ph/ 02020391. 

[7] H. P. Miles, M. Srednicki and D. Wyler, Phys. Lett. B 120 (1983) 346; 
J. P. Derendinger and C. A. Savoy, Nucl. Phys. B 237 (1984) 307; 
M. Drees, Int. J. Mod. Phys. A 4 (1989) 3635; 

J. R. Ellis, J. F. Gunion, H. E. Haber, L. Roszkowski and F. Zwirner, Phys. Rev. D 

39 (1989) 844; 

T. Elliott, S. F. King and P. L. White, Phys. Lett. B 314 (1993) 56 [hep-ph/9305282|; 
Phys. Rev. D 49 (1994) 2435 [hep-ph/9308309|; 

U. Ellwanger, M. Rausch de Traubenberg and C. A. Savoy, Phys. Lett. B 315 (1993) 
331 | |hep-ph/9307322| |; 

P. N. Pandita, Phys. Lett. B 318 (1993) 338; Z. Phys. C 59 (1993) 575; 
B. R. Kim, A. Stephan and S. K. Oh, Phys. Lett. B 336 (1994) 200. 

[8] F. Franke and H. Fraa s, Z. Phys. C 72 ( 1996) 309 [ hep-ph/9511275| |; Int. J. Mod. 

Phys. A 12 (1997) 479 | |hep-ph/ 9512366 1; 

F. Franke, H. Fraas and A. Bartl, Phys. Lett. B 336 (1994) 415 (hepjph/9408217|. 

[9] S. Y Choi, D. J . Miller and P. M. Zerwas, Nucl. Phys. B 711 (2005) 83 
| |hep-ph/ 04072091. 



12 



[10] U. Ellwanger and C. Hugonie, Eur. Phys. J. C 5 (1998) 723 |hep-ph/9712300|; Eur. 
Phys. J. C 13 (2000) 681 |hep-ph/9812427|. 

[11] S. Hesselbach, F. Franke and H. Fraas, Phys. Lett. B 492 (2000) 140 

|hep-ph/ 00073101; 

F Franke and S . Hesselbach, Phys. Lett. B 526 (2002) 370 |hep-ph/0 1112851; 
|hep-ph /0210363l 

[12] G. Moortgat-Pick, S. Hesselbach, F. Fra nke and H. Fraas, |hep-ph/ 9909549} 

S. Hesselbach, F. Franke and H. Fraas, hep-ph/0003272: Eur. Phys. J. C 23 (2002) 
149 [hep-ph/0107080|. 

[13] W. Porod, Comput. Phys. Commun. 153 (2 003) 275 lhep-ph/030llSH ; 
R. Lafaye, T. Plehn and D. Zerwas, hep-ph/ 0404282, see also 

http: //sf itter . web. cern. ch/SFITTER/; 

P. Bechtle, PhD Thesis, University of Hamburg, July 2004; 

P. Bechtle, K. Desch and P. Wienemann, hep-ph/0412012 , see also 

http : / / www-f lc . desy . de/ fittino/. 

[14] K. Desch, J. Kalinowski, G. Moortgat-Pick, M. M. Nojiri and G. Polesello, JHEP 
0402 (2004) 035 [hep-ph/0312069|; 

B. C. Allanach et al. [Beyond the Standard Model Working Group Collaboration], 

hep -ph/0402295) 

G. Moortgat-Pick, hep-ph/0406180 



[15] U. Ellwanger, J. F. Gunion and C. Hugonie, JHEP 0502 (2005) 066 
|hep-ph/ 0406215[. 

[16] P. Skands et al, JHEP 0407 (2004) 036 Mhep-ph/03lTT2"3") . 

[17] L. Randall and R. Sundrum, Nucl. Phys. B 557 (1999) 79 | |hep-th/9810"T55| ; 

G. F. Giudice, M. A. Luty H. Murayama and R. Rattazzi, JHEP 9812 (1998) 027 

|hep-ph /9810442[; " ' 

A. Pomarol and R. Rattazzi, JHEP 9905 (1999) 013 [hep-ph/ 9903448 [. 

[18] A. J. Barr, C. G. Lester, M. A. Parker, B. C. Allanach and P. Richardson, JHEP 0303 
(2003) 045 |hep -ph/02 08214[. 

[19] C. Hensel, DESY-THESIS-2002-047. 

[20] J. F. Gunion and S. Mrenna, Phys. Rev. D 64 (2001) 075002 [he p-ph/0103167| ; 

[21] G. Belanger, F. Boudjema, A. Pukhov and A. Semenov, hep-ph/ 0405253 

[22] U. Ellwanger and C. Hugonie, Eur. Phys. J. C 25 (2002) 297 |hep-ph/ 9909260[; 
D. J. Miller, R. Nevzorov and P. M. Zerwas, Nucl. Phys. B 681 (2004) 3 
piep-ph/0304049[. 



[23] D. J. Miller and S. Moretti, hep-ph/0403137 



13 



[24] U. Ellwanger, J. F. Gunion, C. Hugonie and S. Moretti, hep-ph/ 0305 109: 
J. R Gunion and M. Szleper, hep-ph/ 0409208, 

[25] For details see the Polarization Working Group report of the TOWER group', 

http : / / www . ippp . dur . ac . uk/ ~ gudrid/ power. 

[26] H. E. Haber and G. L. Kane, Phys. Rept. 117 (1985) 75. 

[27] G. Moortgat-Pick, A. Bartl, H. Fraas and W. Majerotto, |Kep-ph/ 002253 . 

[28] The ATLAS Collaboration, ATLAS TDR 15, CERN/LHCC/99-15 (1999); 

M. Chiorboli and A. Tricomi, CMS-NOTE-2004-029; 

B. K. Gjelsten, D. J. Miller and P. Osland, JHEP 0412 (2004) 003 [hep-ph/0410303|; 
G. Weiglein et al. [LHC/LC Study Group], |hep-ph/0410364| 

[29] B. C. Allanach et al, Eur. Phys. J. C 25 (2002) 113 [hep-ph/0202233|; 
N. Ghodbane and H. U. Martyn, |hep-ph/0201233P 
see also http : / / www . cpt . dur . ac . uk/ "georg/ sps/. 

[30] G. Polesello, J. Phys. G 30 (2004) 1185. 

[31] G. J. Gounaris, J. Layssac, P. I. Porfyriadis and F. M. Renard, Phys. Rev. D 71 (2005) 
075012 p\ep-ph/0411366| |. 

[32] H. Fraas, F. Franke, S. Hesselbach, G. Moortgat-Pick and G. Polesello, in prepara- 
tion. 

[33] G. Polesello, private communication. 



14 



